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Abstract: The family symmetry £77(3) <8> 27" (1) is proposed to solve flavor problems 
about fermion masses and flavor mixings. It's breaking is implemented by some 
flavon fields at the high-energy scale. In addition a discrete group Z 2 is introduced 
to generate tiny neutrino masses, which is broken by a real singlet scalar field at 
the middle-energy scale. The low-energy effective theory is elegantly obtained after 
all of super-heavy fermions are integrated out and decoupling. All the fermion 
mass matrices are regularly characterized by four fundamental matrices and thirteen 
parameters. The model can perfectly fit and account for all the current experimental 
data about the fermion masses and flavor mixings, in particular, it finely predicts 
the first generation quark masses and the values of 0[ z and Jq P in neutrino physics. 
All of the results are promising to be tested in the future experiments. 
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I. Introduction 

The precise tests for the electroweak scale physics have established plenty of 
knowledge about the elementary particles pp. The standard model (SM) has been 
evidenced to be indeed a very successful theory at the current energy scale [2]. 
However, there are some imperfections in the SM, among other things, a ugly defect 
is too many parameters exist in the Yukawa sector. As a result, fermion masses 
and flavor mixings seem intricate and ruleless. During the past decade a series 
of new experiment results about B physics and neutrino physics have told us a 
great deal of information about flavor physics [3]. What deserves to be paid special 
attention are some facts as follows. The mass spectrum of quarks and charged 
leptons emerges a large hierarchy, which ranges from one MeV to a hundred GeV 
or so pQ. The neutrinos have been verified to have nonzero but Sub-eV masses 
[I], nevertheless, that their nature is Majorana or Dirac particle has to be further 
identified by experiments such as 0uf3f3 [5]. On the other hand, the flavor mixing in 
the quark sector is distinctly different from one in the lepton sector. The former has 
small mixing angles and its mixing matrix is close to an unit matrix [B] , whereas the 
latter has bi-large mixing angles and its mixing matrix is close to the tri-bimaximal 
mixing pattern [7J. In the lepton mixing, it is yet in suspense whether sinOiz is zero 
and the CP violation vanishes or not [8]. These impressive puzzles always attract 
great attention^, and also are expected to be explained by new theories beyond the 
SM. The issues in the flavor physics possibly implicate great significance. They are 
not only bound up with origin of matter in the universe [10], but also in connection 
with the genesis of the matter-antimatter asymmetry and the original nature of the 
dark matter [TT] . 

Any new theory beyond the SM has to be confronted with the various intractable 
issues mentioned above, however, some approaches and theoretical models have been 
proposed to solve them [12]. For instance, the Froggatt-Nielsen mechanism with 
£7(1) family symmetry can account for mass hierarchy [13]. The discrete family 
group can lead to the tri-bimaximal mixing structure of the lepton mixing ma- 
trix [H]. The non-Abelian continuous group SU(3) is introduced to explain the 
neutrino mixing [T5]. In [IE], the model with the family group SO (3) successfully 
accommodates the whole experimental data of quarks and leptons. In addition, 
some models of grand unification (GUT) based on 5*0(10) symmetry group can 
also give some reasonable interpretations for fermion masses and flavor mixings [TT] . 
Although these theories seem successful in explaining some of the flavor problems, 
it seems very difficult for them to solve the whole flavor problems all together. It 
is especially hard for some models to keep the principle of the smaller number of 
parameters. On all accounts, it now remains to be a large challenge for theoretical 
particle physicists to uncover these mysteries of the flavor physics. 

As remarked above, the flavor puzzle is very difficult and complex to solve com- 
pletely. In this works, we consider a new approach and construct a model with 
fewer parameters to try to understand these problems. On the one hand, we believe 
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that there is some inherence relations among all kinds of fermion mass and mixing 
parameters. The family symmetry SU(3)p, which is a family symmetry group of 
the three generation fermions, is appropriate for seeking the relations. In addition, 
we introduce some super-heavy fermion and flavon fields which appear only at the 
high-energy scale. They communicate with the low-energy fermions of the SM by 
flavor interaction. The different super-heavy fermions and flavons are distinguished 
by an appended Abelian group U(1)n. On the other hand, the model with two 
Higgs doublets and three generation right-handed neutrino singlets is theoretically 
well- motivated extension of the SM[18]. On that basis, we introduce a real scalar 
singlet and append a discrete group Z 2 , under which these scalar and right-handed 
neutrino singlets all reverse sign. The model goes through three steps of breakings. 
Firstly, the family symmetry SU (3)f®U(1)n is broken at the flavon dynamics scale. 
It is accomplished by means of the flavon fields developing the vacuum structures 
along the specific directions. After all the super-heavy fermions are integrated out 
and decoupling, the low-energy effective theory is elegantly obtained. Secondly, the 
discrete symmetry Z 2 is broken at the middle-energy scale by the real scalar devel- 
oping non- vanishing vacuum expectation value (VEV). This directly causes that the 
effective Yukawa couplings of neutrinos are drastically suppressed, thus far smaller 
than ones of the other fermions. This also becomes a source of the neutrino tiny 
masses. Lastly the electroweak symmetry breaking is completed, all of quarks and 
leptons attain Dirac masses. All the fermion mass matrices are regularly given and 
characterized only by the four fundamental matrices and fewer parameters. Finally, 
the model can naturally and correctly give rise to the fermion mass spectrum and 
flavor mixing angles. All the numerical results are very well in agreement with the 
current experimental data. 

The remainder of this paper is organized as follows. In Section II we outline 
the model. In Sec. Ill, the symmetry breaking procedures are introduced and the 
fermion mass matrices are discussed. In Sec. IV, we give the detailed numerical re- 
sults about the fermion masses and flavor mixings. Sec. V is devoted to conclusions. 

II. Model 

The model is based on the symmetry group SU (3)c®SU(2)l®U(1)y ®SU (3).p<g> 
U(1)n ® Z2, among them, the first three subgroups are namely the SM symmetry at 
the low-energy scale. The family symmetry at the high-energy scale is characterized 
by the subgroups SU(3)f ® U(1)n- The subgroup Z 2 is a discrete symmetry at the 
middle-energy scale. The model particle contents and their quantum numbers under 
the family symmetry subgroups are listed in the following. The low-energy fermions 
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and Higgs scalar fields consist of 

Q L ~(3,0), u R ~(3,0), d fl ~(3,0), 

L L ~(3,0), i/ fl ~(3,0), e fl ~(3,0), 

i?i~(l,2), iy 2 ~ (1, -2) , 0~(1,O). (f) 

The three generation of fermions are in 3 representation of the family subgroup 
SU(3)f, and they have no charges of the subgroup U(1)n, and so on. Under the 
SM group, the representations of these fields are clear as usual but only ur and (f> 
are singlets. 

We introduce some super-heavy fermion fields as follows 

^42 
~(3,2), % ' ~ (3, -) , r/ 3 ' ~(3,g), r/ 4 ' ^^'3)' 

%'^(3j), ^"-(3,1), #"~(3,±), C~(M), X"~(3,0), 

r/f e ~(3,-2), r/f~(3,-^), ^-(3,-^), r/f ~(3,-|), C* e ~ (1, -1) , 

(2) 

whose left-handed and right-handed fields are unified. The superscripts respectively 
indicate the corresponding right-handed fermions in (1), namely under the SM group 
the quantum numbers of the super-heavy fermions are the same as ones of the low- 
energy right-handed fermions. In comparison with the super-heavy quark fields, 
the super-heavy lepton fields have \ u instead of ( u . These super-heavy quarks and 
leptons are possessed of the super-heavy masses, so they appear only in the very 
high energy circumstances. 

We also introduce the super-heavy scalar flavon fields such as 

Fi ~ (8, i) , F 2 ~ (8, |) , F 3 ~ (8, jj) , F 4 ~ (8, l -) , 

7\ ~ (3, 1) , T 2 ~ (3, -1) . (3) 

All of them are singlets under the SM group, but under the family subgroup, 
Fi, ■ ■ ■ , F4 are hermitian octet representations, and T\, T2 are complex triplet rep- 
resentations. In addition, they have different charges of U(1)n- These flavon fields 
are responsible for the family symmetry breaking. 

Finally, we define the discrete group Z 2 as follows. Only the v R and <fi fields are 
transformed as 

Ur > -U R ,0 > -<j> , (4) 

and all of the other fields are uniformly transformed as themselves. 
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Under the model symmetry group, the gauge invariant Yukawa couplings in the 
quark sector are written as 

& q = QIH2V1R + vTl (TiQ + F x rf 2R + Fsrfa + F^ R ) + ^f L F,rfi R 

+ WlFA + WlFa^ r + {QTl + ljf L F 2 + 7^F 3 + ^F 4 ) u R 

+ QlH lV d m + < (T 2 ( d R + Fl4 R + F; V d m ) + WlF;4r 

+ (Cffi + WlF; + WlF;) d R + h.c. . (5) 

For the sake of concision, we have left out the coupling coefficient at the front of 
each term in (5), which should be ~ Easy to notice, the couplings in the 

up-type sector are different from ones in the down-type sector. Those F-type flavon 
fields in the down-type sector are complex conjugate form. The Yukawa couplings 
in the lepton sector are similarly given as 

ft = LlH 2 rT 1R + Wl (FiVm + % + F A r£ R ) + ^F 3V » 4R 

+ WlFavIr + WlFaVir + {vkF 2 + 7j^F 3 + 7fcF 4 ) X r + Yl^r 

+ T L E x rf 1R + Wl (T2CR + F? V e 2R + F; V e 3R ) + WlF^ 

+ (ClTl + WlF* 2 + WlF;) e R + h.c. , (6) 

likewise, all the coupling coefficients are omitted. In comparison with the quark 
sector, the lepton sector has the exclusive terms related to x u instead of These 
differences play key roles in generating distinct masses and mixings for quarks and 
leptons. The (5) and (6) Lagrangian indicate that the flavor interactions among 
the super-heavy fermions are transmitted by means of the super-heavy flavon fields, 
and the effects is ultimately transferred to the low-energy fermions after the multiple 
transmissions. 



III. Symmetry Breakings and Fermion Mass Matrices 

The model symmetry breakings go through three stages. The first step of the 
breaking chain is that the subgroups SU(3) f®U(1)n break to nothing, namely the 
family symmetry vanishes. This is implemented by the flavon fields iq, • • • , F 4 , T 1; T 2 
developing VEVs along specified directions in the family space. The detailed vacuum 
structures are as follows 

<*i> 
Af 

(31 

(7) 




(3A 3 - v^Ae + v^Ag) , 
/FA 
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where A F is the family symmetry breaking scale, that is the dynamics scale of the 
super- heavy fermions and flavons, which is usually close to Planck scale of 10 19 GeV. 
The only breaking parameter e is a ratio of the F-type VEV to the T-type VEV. 
We consider that the former is one order of magnitude smaller than the latter, thus 
e is ~ ^(0.1). The matrices Ai,A2,--- , Xs are the standard Gell-Mann matrices 
representing the generators of 577(3). As before one coefficient of ^(1) is implied 
in the right formula of each wave notation in (7). Below the scale , all of the 
flavon fields develop the vacuum states with the structures, consequently the family 
symmetry is broken. It can be seen from (7) that the breakings of 7\ and T 2 bring 
the family symmetry down from SU(3) to SU(2). On the other hand, the breakings 
of Fi and F 2 occur along the direction of the subgroup 52(2 -h- 3) in the family space, 
which is a permutation group between the second generation fermions and the third 
ones, and the F 3 and F 4 breakings are respectively in the directions of the subgroup 
^(l -H- 3) and S2 (1 -H- 2). The multiple breakings lead the family symmetry to 
disappear completely. This set of breaking patterns can be determined in principle 
by the self-interaction potential of every flavon field, but we do not go into the 
details. We directly adopt the intuitive and reasonable scenarios, in particular, it 
turns out to be a great success in fitting experimental data. 

When the energy scale goes down to the value far smaller than Ap, all the super- 
heavy fermions are actually decoupling. After all of them are integrated out from 
the original Lagrangian, then an effective Yukawa Lagrangian at the low energy is 
derived as 

^Yukawa = Q~lH 2 Y u u r + Ql^Y^R + T L E x Y e e R + T7 L H 2 Y u ^v R + h.c. (8) 
with Yukawa coupling matrices 

Yu = ViRi + Vl e 2 R2 + 2/ 3 M ^3 + Vt e A R, , 
Y d ^y d 1 R 1 + yls 2 R 2 + y d 3 e 3 R* 3 , 
Y e = y e 1 R 1 + y e 2 e 2 R 2 + y e 3 e 3 R* 3 , 

Y v = y\ e 2 R 2 + y v 2 e 3 R 3 + y v 3 e 4 R 4 , (9) 

where 














Ri = 


(: 





I) 



















( 





1-i 


Rz = 


V2 \ 




+ i 



1 + 
1-i 






' 


-1 


-V2 


-V2 








R4 


.(i 




\ 



1 I . (10) 


In (9), y", • • • ,y 3 are some effective coupling coefficients, which are left out before, 
now they are visibly retrieved and written out. These coefficients are mostly ~ ^(1), 
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Figure. 1. The graph of generating the effective Yukawa coupling matrices from the 
family symmetry breaking, (a) and (b) respectively give rise to the terms of Ri and 
R 2 in Y u . 

we take them as real numbers without loss of generality, so all the Yukawa matrices 
are hermitian. We can illustrate this procedure, for instance, the terms of Ri and 
R 2 in the Y u matrix are generated respectively by (a) and (b) in the Figure 1, 
and so on. The effective theory, which is valid until the scale A^, includes the 
SM fermions as well as three singlet right-handed neutrinos, moreover, has two 
doublet and one singlet Higgs fields. Since the term involving vr and <fi in (8) is 
drastically suppressed by Ap, it is far smaller than the other terms in Lagrangian, 
furthermore, it is also non-renormalizable. It is very clear from (9) and (10) that this 
set of Yukawa coupling matrices indeed have some regular and intrinsic relations. 
Several notable characteristics can be seen very easy. First, every Yukawa matrix 
is a linear combination from the four fundamental matrices Ri, ■ ■ ■ ,R^, and the 
combination coefficients are expanded by a power series of e. By virtue of such 
structures the elements of every Yukawa matrix show themselves large hierarchy. As 
a result, the R\ and R 2 terms, namely the e° and e 2 terms, will respectively dominate 
the third and second generation fermion masses. The rest of the terms will make 
main contributions to the first generation fermion mass. Second, in contrast with 
Y u ,Y d ,Y e , the leading term of R\ is no in Y u . In view of these structure features of 
the Yukawa matrices, it is in the course of nature that the transformation matrices 
diagonalizing Y u , Yd, Y e are all close to the unit matrix, whereas the transformation 
matrix diagonalizing Y v is approximately the tri-bimaximal mixing pattern. This 
difference is the principal source of generating distinct flavor mixings for the quarks 
and leptons. Third, the imaginary elements of R3 are the only source of the C 
and CP violations in the Yukawa sector. To sum up, the four i?-type matrices and 
the e parameter all together make up the skeleton frame of every Yukawa matrix, 
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therefore they play key roles in fermion masses and flavor mixings. 

The second step of the breaking chain is that the discrete subgroup Z 2 is broken 
by the real singlet scalar field (p developing VEV as follows 

A^ = K - 

The breaking scale is considered as some intermediate value between the family 
breaking scale and the electroweak breaking scale. If the parameter k is about 
10 -10 or so, the effective Yukawa couplings of the neutrinos will be drastically sup- 
pressed owing to the k factor, thus they will be far smaller than the Yukawa couplings 
of the charged fermions. This possibly becomes a source of the neutrino tiny masses, 
of course, it is different from the usual see-saw mechanism [T9] . 

After the Z 2 breaking, the model remaining symmetry is exactly the SM symme- 
try group. The last step of the breaking chain is namely the electroweak symmetry 
breaking. It is accomplished by the doublet Higgs fields Hi and H 2 developing VEVs 
as follows 



(Hi) ( \ (H 2 ) ( sin/3 



cos/3 J' Vew V 



(12) 



where v ew is the electroweak scale and tan/3 is a ratio of the up-type VEV to the 
down-type VEV. After the electroweak breaking all of the quarks and leptons obtain 
Dirac masses. The whole fermion mass terms are now given as 

- -^mass = uZM u u R + d L M d d R + eZM e e R + 17[M u v R + h.c. (13) 

with the mass matrices 

M u = -v ew sin/3 Y u , M v = -« v ew sin/3 Y v , 

M d = -v ew cos(3 Y d , M e = -v ew cos(3 Y e . (14) 

Now three new parameters K,v ew ,tan(3 are added into the model besides the e 
and y-type parameters in (9). All of these quantities are undetermined except the 
electroweak scale v ew . However, some parameters among them are in the form 
of product factors in the mass matrices, there are actually three non-independent 
parameters. We have some freedoms to remove the non-independent parameters, for 
example, this three parameters , yi, Vi can be absorbed collectively by redefinitions 
of the three parameters e, tan/3, k, so each of them will be equal to one instead of 
free parameters hereinafter. Therefore, the effective independent parameters in the 
model are only thirteen in all. It can be seen from (14) that v ew dominates mass 
scale of the quarks and the charged leptons, tan/3 is responsible for mass split of the 
up-type and down-type fermions, and the factor k causes that the neutrino masses 
are far smaller than ones of the charged fermions. Since this three parameters are 
only some product factors in the mass matrices, anyway, they have no influence on 
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the flavor mixings. The mass hierarchy and the flavor mixings are still controlled 
mainly by the s parameter and the f?-type matrices. In a word, this set of mass 
matrices properly embody all of the information about fermion mass hierarchy, flavor 
mixings and the CP violations. 

In virtue of the model's intrinsic characteristics, all the fermion mass matrices 
are hermitian, therefore all of fermion mass eigenvalues are conveniently solved by 
diagonalizing them as follows 

Ul M u U u = diag (m u , m c , m t ) , U ] v M v U u = diag (mi,m 2) m 3 ) , 

U\ M d U d = diag (m d , m s , m b ) , U\ M e U e = diag (m e , m M , m T ) . (15) 

Because the R\ and R2 terms are respectively the leading and next-to-leading terms 
in the mass matrices, the second and third generation of the quark and charged 
lepton masses can be calculated approximately such as 

m c « v ew sin/3 (e 2 + e 4 (2 + yf - yl)) , m t « v ew sin(3 (yl - e 2 ) , 

m s w v ew cos/3 (^e 2 + e 4 ^ + yf J j , m b w v ew cosj3 (-yf + e 2 ) , 

m M » v ew cos/3 [s 2 y e 2 + e 4 + ^) ^ , m r w v ew cos(3 (yl - e 2 y c 2 ) . (16) 

However, the first generation of the quark and charged lepton masses have no such 
approximate expressions since they depend on all the terms of every mass matrix. 
It can be seen from (16) that, in the leading approximation, there are the mass 
relations 

m c e 2 m s e 2 m a e 2 , . 

m t yf — £ nif, —y? + e A m T sx _ £ 2 

We can easily estimate values of some parameters from (16) and (17). Finally, the 
flavor mixing matrices for the quarks and leptons are respectively given by [20J 

Ul U d = Uckm , Ul U v = U PMNS . (18) 

The mixing angles and CP- violating phases in the two unitary matrices of Uckm 
and Upmns can be worked out by the standard parameterization in particle data 
group p. 



IV. Numerical Results 

Now we present the model numerical results. As is noted earlier, altogether the 
model parameters involve the ten y-type coefficients, the three breaking parameters 
£,K,tanf3, and the electroweak scale v ew . Once this set of parameters are chosen 



9 



as the input values, according to the model we can calculate the various output 
values of the fermion masses and flavor mixings, moreover, all of the results can be 
compared with the current and future experimental data. 

The electroweak scale v ew is essentially determined in the gauge sector by weak 
gauge boson masses and gauge coupling constant. The accurate measures have 
given v ew = 174 GeV. The other thirteen parameters are really free parameters in 
the Yukawa sector, however, all of them have to be fixed by fitting the experimental 
data of the fermion masses and flavor mixings. Because the number of the model 
parameters is much less than the experimental values of the masses and mixings, 
and the majority of them have precisely been measured, the space of the model 
parameters is constrained very narrow and the tuning scope of the parameters is 
indeed very small. Although the fit is a non-trivial and no easy one, in advance one 
can find some parameter values by (16) and (17), and then the global fit can be 
successfully finished. We here give the values of the best fit instead of the detailed 
numerical analysis. The input values of the model parameters are elaborately chosen 
as follows 

£ = 0.0845, k = 1.83 x 10~ u , tan(3 = 12.37, 

y? = l, $ = -1, y% = -0.7, = -0.292, ^ = 2.58, 

j£ = 3.62, j& = 7.3, ^ = 0.133, y 2 e = 0.96, y% = 0.871 . (19) 

These values are reasonable and consistent with the previous estimates. Each of the 
y-type parameters is dedicated to a certain impact on the fermion masses and flavor 
mixings, for instance, y% has main impact on the quark sin9f 3 and CP- violating 
phase, while yf makes main contribution to the quark sin9\ 2 , and so forth. 

Finally, a variety of the numerical results calculated by the model are in detail 
listed in the following. For the quark sector, all of mass eigenvalues and mixing 
angles are (mass in GeV unit) 

m u = 0.00246 , m c = 1.272 , m t = 172.2 , 
m d = 0.00472, m s = 0.101 , m b = 4.2 , 

sf 2 = 0.2255 , s| 3 = 0.0414 , sf 3 = 0.00342 , 6 q = 0.375 rr w 67.5° , (20) 

where s a p = sin8 a/ s, in addition, the Jarlskog invariant measuring the CP violation 
is calculated to 

J£ P « 2.87 x 10~ 5 . (21) 

It is very clear that the above results are very well in agreement with the current 
measures of the quark masses, mixing and CP violation pQ. Although the first 
generation of the quark masses have not been accurately measured so far, their 
values are finely predicted to be about the center values of the experimental limits. 
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For the lepton sector, the parallel results are 

m e = 0.511 MeV, m M = 105.7 MeV , m T = 1778MeV, 

mi = 0.172 x 10" 2 eV, m 2 = 0.888 x 10~ 2 eV , m 3 = 5.01 x 10~ 2 eV , 

s[ 2 = 0.565, 43 = 0.751, s[ 3 = 0.109, 5 l = -0.82 vr « -147.6° . (22) 

The charged lepton masses are completely identical with ones in the particle list [TJ. 
For the sake of comparison with the experimental data, the common used quantities 
in neutrino physics are explicitly calculated as follows 

Am 2 x « 7.60 x 10~ 5 eV 2 , Am 2 2 « 2.43 x 10~ 3 eV 2 , 

sm 2 2^ 2 « 0.869 , sm 2 2^ 3 « 0.984 , sm 2 2^ 3 « 0.047 , 

J^p ~ -0.01334 , (23) 

where Am^ = m 2 — m 2 . These results are excellently in agreement with the recent 
neutrino oscillation data [21] . In particular, the model predicts that the heaviest 
one of the three generation neutrinos is about 0.05 eV, the lepton mixing angle #{3 is 
~ 3.3° but nonzero, in addition, the CP-violating effect is of the order of 10 -2 in the 
lepton sector, which is three order of magnitude larger than one in the quark sector. 
Since all the neutrinos in the model are Dirac-type rather than Majorana-type, the 
neutrinoless double beta decay is inevitably nought. Although these quantities have 
not strictly measured by now, some running and coming neutrino experiments are 
on the way toward these goals [22J . We have confidence that all the predictions are 
promising to be tested in the near future. 

To sum up the above numerical results, in fact, only with the thirteen parameters 
does the model accurately and excellently fit the total twenty values of the fermion 
masses and flavor mixings. All the current measured values are exactly reproduced, 
meanwhile, all the non-detected values are finely predicted in the experimental lim- 
its. All of the results are naturally produced without any fine tuning. This fully 
show a strong prediction power of this model. In the case of the best fit, the sizes 
of the parameters and y% both are coincidently one, the reason about it is yet 
unknown and expected to research deeply. 

V. Conclusions 

In the paper, we have suggested a new model to solve the fermion masses and 
flavor mixings, which is based on the family symmetry SU(3)f <8> U(1)n and the 
discrete group Z 2 . The family symmetry breaking is carried out by means of the in- 
troduced super-heavy fermion and flavon fields. After all of the super-heavy fermions 
are integrated out and decoupling, the low-energy effective theory is obtained with 
the regular Yukawa coupling matrices. The Eq parameter and the four fundamental 
R-type matrices all together make up the skeleton frame of the Yukawa matrices. In 
fact they play leading roles in the model, namely they dominate the fermion mass 
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hierarchy and flavor mixing results. The discrete group Z 2 is broken by the singlet 
scalar field at the middle-energy scale. This leads that the Yukawa couplings of the 
neutrinos are drastically suppressed, and then gives rise to the tiny nature of the 
neutrino masses. That set of the fermion mass matrices derived from the model 
symmetries and their breakings are characterized only by the thirteen effective pa- 
rameters. The model successfully and perfectly fits all the current experimental data 
about the fermion masses and flavor mixings, in particular, it finely predicts the first 
generation quark masses and the values of 0' 3 , Jq P in neutrino physics. All of the 
results are excellent and inspiring, and also fully show a great prediction power of 
the model. Finally, we expect all the results to be tested in future experiments on 
the ground and in the sky. These experiments will undoubtedly provide us more 
important information about the flavor physics, and then enlighten us to understand 
finely the mystery of the universe. 
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